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Abstract: We report on the mesoscale fabrication and characterization of 
polymeric templates for isotropic photonic materials derived from hyper- 
uniform point patterns using direct laser writing in a polymer photoresist. 
We study experimentally the microscopic structure by electron microscopy 
and small angle light scattering. Reducing the refractive index mismatch by 
liquid infiltration we find good agreement between the scattering data and 
numerical calculations based on a discrete dipole approximation. Our work 
demonstrates the feasibility of fabricating such random designer materials 
on technologically relevant length scales. 
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1. Introduction 

Coherent scattering can strongly influence the macroscopic transport of electrons or photons. 
In the limit of weak scattering from periodic lattices the classical concept of Bloch allows to 
distinguish propagating and non-propagating wave vectors k. For a periodic one-dimensional 
structure with a lattice constant d forbidden wave vectors appear around the edge of the Bril- 
louin zone k = n/a, the latter being a signature of back-reflection from Bragg planes with 
scattering vectors q = 2k = 2n / a [1 1. Despite the beauty of Blochs theory it is well known that 
the model covers only a small part of actual electron transport phenomena in metals or semi- 
conductors. A prime example is the fact that many alloys and even liquids display metallic or 
semiconducting properties ]2). 

The broader concept of electronic scattering from ordered or disordered structures can be 
readily transferred to the case of electromagnetic waves as shown in the pioneering work of 
Yablonovitch and John (3]|4]. Following their discovery there have been enormous research 
efforts aimed at manufacturing photonic semi-conducting materials. Again these efforts have 
largely concentrated on crystalline structures, owing to the simple deterministic design rules 
and the fact that band-structure calculations are well established for crystals fl] [5] |6). The 
structural anisotropy of periodic structures, both in real space and reciprocal-space, however, 
imposes physical limitations for three-dimensional photonic materials. A direct consequence 
of anisotropy is the appearance of stop bands with varying strengths depending on the direc- 
tion of wave propagation. In order to achieve complete band-gaps in three dimensions the gaps 
in all directions have to overlap, a condition that is difficult to achieve or that requires a very 
large dielectric contrast. As emphasized previously [7 8], the availability of isotropic pho- 
tonic materials would thus be highly desirable. A step in this direction has been made with the 
discovery that quasi-crystalline structures equally possess photonic properties while having a 



much higher degree of rotational symmetry lr9l [TTl . Quasicrystalline structures are aperiodic but 
nevertheless exhibit sharp Bragg peaks [ 10]. Since they are ordered and can be projected from 
higher-dimensional crystalline structures they can be considered a case in between crystalline 
and amorphous. Disordered dielectric heterostructures have also been considered candidates 
for photonic band gap materials lfT2l [T3l . However, until recently, the design rules needed to 
derive such structures have remained obscure. A number of studies have looked into materials 
with short-range order such as photonic liquids and glasses derived from self-assembled col- 
loidal lfl4l [151 [TBI or biomimetic materials 11171 . Although Bragg-like scattering and pseudo 
gaps have been observed no conceptual proof for the existence of a complete band gap could be 
derived even at arbitrarily high dielectric contrast. Moreover, the interplay between short range- 
order and Anderson localization of light in three dimension, highlighted in the early work of 
Sajeev John [4|, remains an unsolved question. 

Recently Florescu, Torquato and Steinhardt have taken a different route to design amorphous 
structures with a full photonic bandgap (PBG)||7]|. The authors make the general claim that 
three structural conditions have to be satisfied to obtain a PBG: hyperuniformity, uniformlo- 
cal topology, and short-range geometric order. Structures that fulfill these conditions display a 
vanishing structure factor S(q) at small but finite wavenumbers |q| < |q c |, a structural property 
coined stealthy hyperuniformity. Above a threshold dielectric contrast such structures are pre- 
dicted to display a full PBG. A trivial finding is that all crystalline and quasicrystalline materials 
fall into this class. Their main point however is that they could identify peculiar structures that 
are disordered and thus completely isotropic but still fulfill the above mentioned conditions Q. 
It has been shown that for exampled randomly jammed packings of spheres posses such hype- 
runiform long-range correlations 1181 . In addition to the underlying geometric order, the local 
topology plays an important role HJ. Mapping hyperuniform point patterns with short-range 
geometric order into tessellations allows the design of interconnected networks that gives rise 
to enhanced photonic properties. Numerical calculations for such two- and three-dimensional 
disordered designer materials [7 , 8 1 indicate the presence of a robust full band-gap in the limit 
of sufficiently high dielectric contrast, typically n > 3 in air. Preliminary experimental data 
obtained from 2D hyperuniform structures in the microwave regime seem to confirm these pre- 
dictions ifnn . 

Here we report on the first experimental realization of micron scale three-dimensional poly- 
meric templates, n ~ 1.52, based on the design rules derived by Florescu and coworkers Q. 
Our present aim is to demonstrate the feasibility of fabricating such random designer materi- 
als on technologically relevant length scales, a first major step towards the fabrication of three 
dimensional amorphous full PBG materials. 

2. Results 

2.1. Nanofabrication of designer disordered materials in a polymer photoresist. 

We fabricate mesoscale polymer structures by direct laser writing (DLW) (Photonic Profes- 
sional, Nanoscribe, Germany). The method allows to replicate a 3D structure voxel by voxel 
into a polymer photoresist with submicron resolution using a two-photon polymerization pro- 
cess [6|. We first implement the protocol suggested by Florescu and coworkers to map a hype- 
runiform point pattern into tessellations for photonic materials design |7|. As a seed structure 
we use the centroid positions from a maximally randomly jammed assembly of spheres of di- 
ameter a = 3.31/im with a volume filling fraction of (j) ~ 0.64 [20 1. As shown by Torquato 
and coworkers [18] such structures indeed possess the required hyperuniform long-range cor- 
relations. Locally the seed pattern displays pronounced short-range correlations due to close 
packing of the original assembly of spheres. Therefore the real space pair correlation func- 
tion g(r) is sharply peaked at contact |r| = d and thus the structure factor S(q) is peaked at 



Fig. 1. Design of amorphous photonic structures based upon a tetrahedral network of el- 
liptical rods, (a) Enlarged view of part of the network obtained by tessellation of a hy- 
peruniform seed pattern, rod cross section 840 x 280nm 2 . (b) close-up showing also the 
grid points lying inside or at the surface of the rods used for the numerical analysis of the 
scattering pattern. 

I lima* = 27t/d ||2TI . The latter is closely related to the well-known diffraction ring observed in 
experiments on dense molecular and complex liquids or glasses ll22ll . Next we perform a 3D 
Delaunay tessellation of the spheres center positions. In this scheme tetrahedrons are formed in 
such a way that no sphere center is contained in the circumsphere of any tetrahedron in the tes- 
sellation. The center-of-mass of neighboring tetrahedrons are then connected resulting in a 3D 
random tetrahedral network with the desired hyperuniform properties. The only parameters not 
yet fixed are related to the shape and cross section of the dielectric rods replacing the fictitious 
connection lines (Fig.[T](a)). 

We find optimal conditions in the writing process when the rods are written into IPG-780 neg- 
ative tone photoresist (Nanoscribe, Germany). The laser writing pen has an elliptical cross sec- 
tion of about 840 x 280 nm 2 . When writing thinner rods the structure becomes mechanically 
unstable while thicker rods lead to overfilling of the structure. The aspect ratio of ca. 3.0 of the 
rods is dictated by the point spread function of the illuminating microscope objective during 
the writing process which is set by the refractive index of about 1 .52 of the material and the 
numerical aperture of 1 .4. For the parameters chosen we estimate the volume-filling fraction of 
the rods to be roughly 8 %. 

In comparison to the well-established case of periodic rod assemblies @ the fabrication of 
these designer disordered structures is very demanding and required several months of opti- 
mization in order to obtain the results reported here. Particular care must be taken how to set up 
the writing protocol and moreover all writing parameters have to be optimized in order to create 
mechanically stable structures. We note that precision direct laser writing is almost always oper- 
ated by sequential writing of lines in 3D (with 840 x 280 nm 2 cross section in our case) and not 
by a layer-by-layer process, commonly employed in lower-resolution commercial 3D printing 
techniques used for rapid prototyping or manufacturing. Therefore, in our fabrication process 
particular care must be taken how to set up the writing protocol since there exist no obvious 
rules how to write a random free-standing network structure at optimal resolution. Moreover it 
must be ensured that the structure remains mechanically stable in a soft gel photoresist through- 



out the writing process of about 1 h. This task is further complicated by intrinsic mechanical 
stresses created upon exposure of the photoresist. The latter leads to substantial deformations 
if the written lines are not attached within seconds to a mechanically rigid superstructure. To 
overcome these challenging problems, encountered in our initial fabrication attempts, we have 
developed a optimized writing protocol. We first divide the entire volume in cubic sub-volumes 
of side length roughly 1 .5D (where D is the average distance between nearest points in the un- 
derlying point pattern). We sort the rods in such a way that i) once all rods belonging to a certain 
cube are written, we proceed to a neighboring cube, ii) first all cubes closest to the substrate 
are filled with rods, then the ones lying higher above and so on until the whole network has 
been written. Fig. [2] displays a representative set of electron micrographs of fabricated struc- 
tures. We succeeded to write structures with either a square (65 x 65/im 2 ) or circular (diameter 
65 jim) footprint with heights varying between h = 4 — 12/im. Structures higher than 4/zm were 
surrounded by a massive wall (Fig.|2]l for enhanced mechanical stability |f6l . 




Fig. 2. Electron micrographs of fabricated hyperuniform three-dimensional disordered 
structures, (a) Normal view of a structure with height h = 8/xm and inner diameter 
d = 65 \lm (b) close-up view (c) focused ion beam cut of the same structure. 



2.2. Structural analysis by light scattering 

We characterize the properties of our samples by measurements of the scattering patterns using 
visible light. For the optical characterization of our sub-mm sized samples we have built a 
small angle scattering instrument consisting of a helium neon laser (A = 632.8 nm), a focusing 
lens (focal length / = 50 mm), two diaphragms to suppress stray light and a white screen, 
positioned at a distance of z = 125mm from the sample, with a central absorbing beam block. 
The scattered light pattern is photographed off the white screen using a digital camera. The 
setup has been calibrated using a small pinhole. Histogram normalization has been applied to 
all images displayed in Fig. [3] This procedure represents a linear transformation of an image 
where the value Pj„ of each pixel is scaled according to [ 24 ] : P out = 255 (P,„ — c)/(d — c). Here c 
and d are the jt-th and ( 100— x)th percentile in the histogram of pixel intensity values. All values 
P out smaller than zero are set to 0, while those larger than 255 are set to 255. In our case x = 0. 1. 
Before calculating radial averages of experimental diffraction patterns, several data processing 
steps were performed. Despite the use of diaphragms we could still observe some stray light 
contributing to the image. This contribution has to be subtracted from the raw data. To this end 
we have acquired an image from a bare glass substrate (empty cell) inserted into the laser beam. 
This image is subsequently subtracted from the raw data. The modulus of the scattering vector 
q is calculated from the radial distance u from the center via the relation tan(0) = u/z and q = 
(4/A) sin(0 /2). These relations apply also for the (toluene or toluene/chlorobenzene) infiltrated 
structures since for small angles the reduction of the wavelength X/n within the sample is 
to a good approximation offset by refraction at the flat sample-air interface. Another small 
correction results from the actual detector acceptance angle being different for each scattering 
angle 0. This gives rise to a correction |23| of the measured data by a factor cos(0)~ 3 . 
Examples for light scattering data recorded for h = 4/im structures are shown in Fig. [3] The 
data clearly reveals a concentric ring profile without Bragg peaks. The pronounced maximum 
indicates short-range order while the ring-shape is a signature of structural isotropy without any 
long-range order. Another feature related to disorder is the speckled appearance of the scattering 
pattern, reminiscent of the laser speckle observed for an arbitrary random structure. Although 
our first observations confirm the overall picture, the ring pattern obtained in air, Fig. [3|a), 
appears blurred with a weaker than expected maximum and substantial low angle scattering. 
Interestingly similar observations have been made for quasicristalline structures ll26ll . In our 
case we can clearly attribute the blurring to multiple scattering. From the attenuation of the 
direct laser beam intensity we can estimate a scattering length of l s ~ 3/im, short even compared 
to the lowest sample studied with h = 4/im. In order to reduce scattering we infiltrate the 
sample first with isopropanol (n = 1.377) and then with toluol (n = 1.496), which leads to a 
gradual reduction of multiple scattering and a sharpening of the diffraction ring, Fig.[3jb),(c), 
while at the same time the ring position remains unchanged. For the latter case the refractive 
index of the polymeric structure is almost matched and the direct beam is attenuated only by 
a few percent signaling the absence of multiple scattering. Similar results are obtained for the 
h = 8/xm structures when using a 1:2 mixture of toluene/chlorobenzene (« = 1.517) as an index 
matching fluid (Fig. |4](b). In both cases a clear and pronounced ring appears in the scattering 
pattern at q ~ 2n/d and for smaller g-values scattering is strongly suppressed. Spatial speckle 
fluctuations can be largely (although not entirely) suppressed by taking radial averages as shown 
in l ig.|4]. 

2.3. Comparison with numerical calculations 

For a quantitative analysis of the experimental results we numerically solve the scattering prob- 
lem using a discrete dipole approximation (DDA) in the single scattering limit [25 1. Given 
the weak scattering contrast in the toluene- (or toluene/chlorobenzene) infiltrated structures we 
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Fig. 3. Small angle light scattering pattern using a X = 632. 8nm laser. (a),(b),(c), scattering 
pattern for a structure of thickness h = 4/im in air, infiltrated with isopropanol (n = 1.377) 
and with toluene (n = 1.496). (d) calculated scattering pattern in the single scattering limit 
for the same structure based on a discrete dipole approximation (DDA). 



expect this approximation to hold very well. First, rods coordinates generated in the manner de- 
scribed above are used to create a 3D binary representation of the network (Fig.[T](a)). Namely 
a sufficiently dense grid is defined and then each grid point is set to 1 if the point belongs to 
a rod and otherwise (Fig. [T] (b)). Using this procedure the lithographic voxel size and el- 
lipsoidal shape (as obtained in our direct laser writing system) can be fully taken into account. 
Consequently, the binary network represents a good approximation of the true structure as man- 
ufactured by laser nanolithography. Following that, a 3D Fast Fourier Transform of this data is 
calculated. The squared modulus of FFT is proportional to the intensity /(q) scattered by the 
structure for a scattering wave vector q = k ko, where k denotes the scattered wave vector 
and ko the incident wave vector. Lines plotted in Fig.[4]have been obtained by averaging over 
several realizations of a structure of a given height h. As shown in Fig.[3](d) the numerical re- 
sults reproduce both the ring diffraction as well as the superimposed random specular structure. 
A comparison of the radially averaged numerical data with experiment, Fig. [4], reveals a very 
good match with no adjustable parameters except for the absolute scale of intensities. 
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Fig. 4. Radially averaged scattering intensity. Symbols: experimental data I(q) for h = 
4/lm, immersed in toluene (n = 1.496), and h = 8/im, immersed in toluene/chlorobenzene 
1:2 (n = 1.517). Solid lines: theoretical calculations for the same heights. 



3. Discussion 

Our optical characterization confirms the high quality of our fabricated samples. Moreover 
the results readily show that due to the finite size the optical properties of the material are 
not yet fully developed. As can be seen in Fig. |4]the peak height increases from h = 4/im to 
h = 8jlim. Preliminary numerical results (data not shown) indicate that the peak height satu- 
rates for heights larger than 20/imm. Moreover both the experimental and numerical results 
show that the extrapolated values for S(q — y 0) are finite. Due to experimental difficulties ac- 
cessing wavenumbers close to the primary beam we are unable to clearly distinguish finite-size 
effects from residual contaminations and experimental artifacts. Similarly the accuracy of our 
numerical calculations for small wavenumbers is limited due to finite size effects. Nevertheless, 
since the seed structure is hyperuniform and since the quality of our polymeric templates is very 
high, we do believe that the polymeric template possesses the contemplated structural proper- 
ties and thus should likely give rise to a full PBG when transferred into high-index dielectric. 
We now turn our attention to the spectral properties of the interconnected network structures. 
Numerical results predict a broad isotropic bandgap in three dimensions for a refractive index 
of n = 3.6 and a volume filling fraction of about 20% [8 |. In the following we discuss the ad- 



ditional processing steps required in order to obtain such strong photonic properties. Building 
upon the successful fabrication reported here the following parameters have to be optimized: 
a) the material refractive index needs to be increased b) the volume filling fraction must be in- 
creased to about 20% and finally the structural features ideally should be further scaled down. 
To realize the first point one can rely on well-established procedures reported in the literature. It 
has been shown that polymeric templates very similar to ours can be transferred into materials 
such as silicon (n = 3.6 for infrared wavelengths) by double inversion retaining the original 
topology [26, 27 1 . Optimal filling fractions have been predicted to be around 20% while cur- 
rently, in our case, the polymer content is only around 8%. This means that for best results one 
either has to increase the polymer volume fraction and perform a double inversion or alterna- 
tively it should be possible to coat directly the polymeric template with a high index material. 
Finally, reducing the structural length scales will require further incremental optimisation of 
the delicate fabrication process. Such experiments are currently underway in our laboratory. 
For the polymeric structures reported here one would expect a mid gap wavelength in air at 
around A ~ 4n/q max ~ 6/im, about four times larger than typical telecommunication wave- 
lengths of X ~ 1.5jttm QQ. Based on preliminary results we expect to be able to reduce the 
typical length scales by at least a factor of two within the next months. We note that if the fea- 
ture sizes are scaled down this also favourably affects the filling fraction as long as the size of 
the laser-writing pen is kept constant. 

4. Summary and conclusion 

We have successfully demonstrated the fabrication of high quality three-dimensional polymeric 
templates for disordered photonic materials. Although a number of further processing steps are 
still required to obtain strong photonic properties, our results already demonstrate the feasibil- 
ity of creating such complex materials on length scales comparable to optical wavelengths. As 
shown previously such polymeric templates can be reliably replicated into materials such as 
silicon. We thus envision a successful transfer of our polymer templates into a high index ma- 
terial within the near future [27|. Moreover, in the present study we have shown that the design 
parameters of polymeric templates can be set precisely using direct laser writing lithography. 
This in turn will allow rigorous experimental testing of the existing theoretical concepts when 
applied to the high index replica. 
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